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INTRODUCTION 

The  goal  of  this  investigation  was  to  determine  the  three-dimensional  structures  of  the  two 
known  human  estrogen  receptors  (ERa  and  ER(3)  complexed  with  receptor-selective  estrogens 
and  antiestrogens  (SERMs).  An  important,  unresolved  issue  is  the  molecular  mechanism(s)  by 
which  SERMs  can  act  selectively  as  full  agonists,  partial  agonists,  or  complete  antagonists  in  the 
control  of  cell  proliferation  and  cell  fate  in  different  tissues.  The  physiological  effects  of  both 
endogenous  and  synthetic  SERMs  are  mediated  by  the  estrogen  receptors  (ERs),  ERa  and  ER{3, 
members  of  the  nuclear  receptor  (NR)  super-family  of  ligand-regulated  transcription  factors. 

The  two  ER  subtypes  have  overlapping  but  distinct  tissue  distribution  patterns  in  vivo  and 
distinct  activation  profiles  at  both  classical  and  complex  promoter  elements.  Thus,  it  is  unclear 
how  ERa  and  ER|3  individually  contribute  to  the  effects  of  known  estrogens.  Ligands  that  act 
differentially  on  ERa  and  ERj3  would  be  potentially  valuable  both  as  tools  to  dissect  the 
biological  roles  of  the  two  ERs  and  as  novel  therapeutics  with  pharmacological  properties 
distinct  from  existing  drugs.  As  part  of  a  search  for  ER  subtype-selective  ligands,  the  synthetic 
compound,  R,R-5,ll-cw-diethyI-5,6,ll,12-tetrahydrochrysene-2,8-diol  (R,R-THC),  was 
identified  as  a  selective  estrogen  agonist  when  bound  to  ERa  and  as  an  antagonist  when  bound 
to  ER(3.  To  better  understand  this  selective  behavior,  a  major  goal  of  this  investigation  was  to 
determine  the  crystallographic  structures  of  human  ERa  and  ER|3  ligand  binding  domains 
(LBDs)  complexed  with  R,R-THC.  These  structures  have  now  been  solved  and  refined, 
suggesting  mechanisms  by  which  this  compound  can  act  as  an  ERa  agonist  and  as  an  ER|3 
antagonist. 

Previously  determined  structures  of  ERa  with  4-hydroxytamoxifen  (OHT)  and  diethylstilbestrol 
(DES)  (1)  revealed  and  defined  a  multipurpose  docking  site  on  ERa  and  ERp  that  can 
accommodate  either  helix  12,  in  the  presence  of  OHT,  or  one  of  several  co-regulators  in  the 
presence  of  DES.  Thus,  agonists  and  antagonists  bind  at  the  same  site  within  the  core  of  the  ER 
LBD  to  induce  distinct  conformations  in  the  transactivation  domain  (AF-2),  especially  in  the 
positioning  of  helix  12.  Interestingly,  R,R-THC  stabilizes  a  conformation  of  the  ERa  LBD  that 
favors  coactivator  association  and  a  conformation  of  the  ER|3  LBD  that  prevents  coactivator 
association.  A  comparison  of  the  two  structures,  combined  with  functional  data,  reveals  that 
R,R-THC  does  not  act  on  ER(3  through  the  same  mechanisms  used  by  other  known  ER 
antagonists.  Instead,  R,R-THC  antagonizes  ERj3  through  a  novel  mechanism  we  term  “passive 
antagonism”.  Paradoxically,  the  R,R-THC-ER(i  structure  is  very  similar  to  the  structure  induced 
by  genistein  (2),  which  acts  as  a  partial  estrogen  through  both  ER  subtypes.  Ongoing 
mutagenesis  studies  have  helped  define  some  of  the  molecular  and  structural  differences  that  are 
responsible  for  these  unanticipated  results.  In  addition,  the  passive  antagonism  mechanism 
suggests  a  novel  approach  to  the  design  of  ligands  that  selectively  antagonize  the  two  ER 
subtypes.  Such  ligands  may  have  novel  therapeutic  properties  that  can  be  exploited  to  prevent  or 
treat  breast  cancer. 

This  investigation  resulted  in  a  major  publication  (3)  (See  Appendix)  that  describes  our  success 
in  achieving  the  goals  stated  in  the  original  proposal.  During  the  third  and  final  year  we 
continued  to  use  site  directed  mutagenesis  to  define  and  characterize  residues  in  the  ERa/(3 
LBDs  that  are  responsible  for  ER  subtype-selective  ligand  recognition.  It  became  clear  that 
ligand  discrimination  by  the  two  ER  subtypes  is  more  complicated  than  first  imagined.  As 
summarized  below,  our  data  illustrate  that  secondary  structure  interactions  and  flexibility  are  key 
determinants  of  ligand  specificity  for  the  two  ER  subtypes,  revealing  novel  aspects  of  nuclear 
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receptor  function.  Equally  important,  residues  outside  the  ligand-binding  pocket  appear  to  be 
critical  for  ligand  discrimination  between  two  closely  related  nuclear  receptors. 

BODY 

The  Statement  of  Work  for  this  grant  listed  four  major  tasks.  These  tasks  are  listed  below  with 
details  of  the  progress  made  toward  achieving  each  goal  during  the  past  three  years. 

Task  1  (First  year):  To  produce  and  purify  bacterially  expressed  ERa-LBD-BJK-1  (R,R- 
THC)  complexes  for  biochemical  analysis  and  crystallization. 

As  described  in  detail  in  the  attached  reprint  (3)  (see  Appendix)  and  in  last  year’s  Annual  Report, 
this  task  has  been  accomplished.  BJK-1  (named  after  Benita  and  John  Katzenellenbogen)  was 
renamed  R,R-THC  once  we  knew  its  chemical  identity. 

Task  2  (Second  year):  Purification,  crystallization,  and  x-ray  analysis  of  ERa/p-LBD-BJK- 
1  and  ERa/p-LBD-ICI  complexes. 

The  major  goal,  which  was  to  solve  the  R,R-THC-ERa/p  LBD  structures,  has  been  achieved,  as 
described  in  the  attached  reprint  (3)  (see  Appendix).  Both  structures  are  discussed  in  detail 
below.  We  did  not  pursue  either  the  ERa-or  ERp-ICI  182,780  complexes  because  the  R,R-THC 
studies  consumed  all  available  personnel  time  for  this  grant.  Also,  the  Hubbard  lab  solved  the 
structure  of  rat  ER{3  LBD  complexed  with  ICI  164,384  in  2001  (4).  It  may  still  be  of  some  value 
to  know  if  the  ERa-ICI  conformation  is  different  from  the  published  ER|3-ICI  conformation. 
Noteworthy  in  the  ICI-ERp  structure  was  the  absence  of  any  association  of  the  transactivation 
helix  (H12)  with  the  LBD,  which  contrasts  with  all  other  known  antagonist-ER  structures,  in 
which  H12  is  associated  with  the  hydrophobic  cleft  formed  by  helices  3, 4  and  5  (Figure  1  and 
Appendix,  Figure  lb).  In  the  ICI-ER(3  complex,  the  terminal  portion  of  the  ICI  bulky  side 
chain  substituent  emerges  from  the  ligand-binding  pocket  and  binds  to  part  of  the  coactivator 
recruitment  site  (AF2)  (Figure  lc),  thereby  preventing  H12  from  adopting  either  its 


Figure  1.  (a)(b)  Orthogonal  views  of  ER|3-LBD  bound  to  the  pure  antagonist  ICI  164,384. 
6*1  Interaction  of  TCT  side  chain  with  coactivator  groove.  Reproduced  from  Pike  et  al.  141. 
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characteristic  agonist  or  antagonist  orientation.  This  complete  destabilization  of  H12  may 
contribute  to  the  more  potent  antiestrogenic  activity  of  both  closely  related  ICI  compounds  when 
compared  to  SERMs  like  tamoxifen  and  raloxifene.  Another  potentially  significant  difference  is 
the  effect  of  the  ICI  164,384  on  the  ER|3  dimerization  interface.  For  all  other  known  ERa/p 
LBD  structures,  the  ligand  has  little  effect  on  the  nature  or  extent  of  the  dimer  interface. 
However,  ICI  appears  to  generate  a  more  open  and  perhaps  less  stable  dimer,  consistent  with 
reports  that  ERa-ICI  dimers  are  less  stable  than  other  ligand-ER  complexes  (5).  It  is  also 
possible  that  a  crystallization  artifact  is  responsible  for  the  observed  ICI-induced  conformation, 
or  that  the  ER|3  dimer  interface  is  more  sensitive  than  ERa  to  ligand  effects  (4).  For  ER(3-R,R- 
THC,  the  dimer  interface  is  very  similar  to  the  ERfTRal  and  ER(3-Gen  structures  (2, 4). 

R,R-THC  ERa/p  structures.  Refinement  of  the  R,R-THC-ERa/p  LBD  structures  continued 
into  the  early  part  of  year  3,  providing  additional  insights  about  key  structural  differences 
between  the  two  ER  subtypes.  The  overall  structures  of  R,R-THC-ERa  LBD  and  R,R-THC- 
ERP  LBD  are  shown  in  Figure  2  and  in  the  Appendix,  Figure  lb.  The  helices  are  arranged  as 
in  Figure  1,  except  that  H12  (shown  in  red)  is  structured  and  either  resting  across  the  LBD 
pocket  (Fig.  2a,  ERa)  or  in  the  hydrophobic  cleft  formed  by  helices  3-5  (Figure  2b,  ERfJ). 
Although  we  were  able  to  obtain  preliminary  crystallographic  solutions  for  both  the  R,R-THC- 


Figure  2.  Overall  structures  of  the  R,R-THC-ER  LBD  complexes.  A.  R,R-THC-ERa 
LBD  +  GRIP1  NR  box  II  peptide.  B.  R,R-THC-ER(3  LBD  (Appendix,  Figure  lb) 

ERa  and  R,R-THC-ERp  LBD  structures  almost  two  years  ago,  ongoing  data  refinement  during 
the  past  year  revealed  a  more  complex  ligand  interaction  with  the  ER|3  pocket  than  proposed 
initially.  As  shown  in  Figure  3,  our  original  assumption  that  fewer  contacts  between  R.R-THC 
and  the  amino  acid  residues  that  line  the  pocket  of  ER(3  were  responsible  for  the  inability  of  this 
complex  to  assume  an  agonist  conformation  was  incorrect.  Instead,  an  equal  number  of  contacts 
occur  between  R,R-THC  and  amino  acids  that  comprise  the  ligand-binding  pocket,  as  compared 
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to  the  R,R-THC-ERa  LBD  complex  (Appendix,  Figure  3b).  The  ERa  complex  structure  has 
now  been  refined  to  a  crystallographic  R  factor  of  20.3%  (Rfree  =  24.2%)  using  data  to  1.95  A 
resolution  and  the  ER|5  complex  structure  has  been  refined  to  a  crystallographic  R  factor  of 
25.9%  (Rfree  =  29.9%)  using  data  to  2.95  A  resolution  (Appendix,  Table  1).  It  proved  especially 
difficult  to  obtain  an  Rfree  value  of  29.9%  for  ERp,  which  is  why  we  did  not  initially  see  the 
additional  contacts  that  were  present  between  R,R-THC  and  the  ER|3  LBD  pocket.  For  example, 
we  did  not  initially  see  the  additional  hydrogen  bond  that  occurs  between  the  A  ring  phenolic  OH 
of  R,R-THC  and  Arg  346,  similar  to  the  interaction  with  Arg  394  in  the  R,R-THC-ERa  structure 
(Figure  3B;  Appendix,  Figure  3b).  Similarly,  additional  R,R-THC  contacts  with  the  ER|3  LBD 
pocket  could  only  be  appreciated  upon  further  refinement  of  the  structural  data.  We  now  feel 
that  16  amino  acids,  rather  than  the  11  first  observed  (Figure  3B;  Appendix,  Figure  3b), 
participate  in  the  stabilization  of  R,R-THC  in  ER0,  which  equals  the  16  amino  acids  that  are 


Figure  3.  R,R-THC  Interactions  with  ER0  LBD  Pocket.  A.  Original  model.  B.  Current 
model  based  on  most  refined  structure  (see  Appendix,  Figure  3b). 

involved  in  the  R,R-THC-ERa  LBD  complex,  although  not  all  residues  are  the  same  (Appendix, 
Figure  3b,  3c).  Also,  a  water  molecule  participates  in  the  hydrogen  bonds  that  stabilize  the  ‘A 
ring’  of  R,R-THC  in  the  ERa  LBD,  but  not  in  the  ER|3  LBD  (Appendix,  Figure  3b,  3c).  The 
R,R-THC-ERa  LBD  pocket  residue  interactions  are  very  similar  to  those  of  the  full  agonists 
estradiol  and  DES  and  nucleate  all  of  the  polar  and  nonpolar  contacts  that  are  required  to  favor 
the  agonist-bound  conformation  of  ERa  in  which  helix  12  rests  across  the  pocket  (Figure  2A; 
Appendix,  Figure  lb). 

Thus,  although  our  hypothesis  of  “passive  antagonism”  has  not  changed,  our  view  of  the 
underlying  mechanism(s)  has  evolved.  As  described  previously,  the  ER|3  LBD  shares  the  same 
overall  fold  as  the  ERa  LBD  (Figure  2;  Appendix,  Figure  lb)  (2).  And  like  the  ERa  LBD,  the 
ERp  LBD  recognizes  R,R-THC  with  a  ligand-binding  pocket  located  within  its  lower  sub 
domain.  However,  the  binding  of  R.R-THC  to  the  ER|3  LBD  does  not  stabilize  the  agonist- 
bound  conformation  of  helix  12  (Figure  2;  Appendix,  Figure  lb).  Also,  helix  12  is  not  bound 
to  the  portion  of  the  coactivator  recognition  groove  formed  by  helices  3, 4,  and  5  as  it  is  in  the 
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OHT-ERct  LBD  (1),  RAL-ERa  LBD  (6),  and  RAL-ERp  LBD  (2)  complexes  (Figure  2; 
Appendix,  Figure  2a),  nor  is  it  disordered  as  it  is  in  the  ICI-ERp  LBD  complex  (Figure  1)  (4). 
Instead,  helix  12  in  the  R,R-THC-ERp  LBD  complex  interacts  with  the  rest  of  the  LBD  in  a 
manner  most  like  that  observed  in  the  GEN-ERp  LBD  complex  (Figure  4;  Appendix,  Figure 
2a)  (2).  Helix  12  is  positioned  similarly,  and  of  comparable  length  (residues  487  to  498  in  the 
R,R-THC  complex  and  residues  487  to  497  in  the  GEN  complex)  in  the  two  structures  (Figure 
4;  Appendix,  Figure  2a). 


Figure  4.  Positioning  of  H12  in  GEN-ERp  LBD  (left)  and  R,R-THC-  ERP  LBD  (right). 

GEN-ERP  LBD  figure  is  reproduced  from  Pike  et  al.  (2). 

Based  on  previous  correlations  between  the  positioning  of  H12  and  the  ability  of  ERa  and  other 
NRs  to  form  coactivator  complexes,  the  conformation  of  the  ERP  LBD  observed  in  both  the 
R,R-THC  and  GEN  complexes  should  be  incapable  of  interacting  with  coactivators  and,  hence, 
should  be  transcriptionally  silent  for  two  reasons.  First,  key  residues  from  helix  12  (Leu  490, 
Glu  493,  and  Met  494)  that  are  predicted  to  form  part  of  the  coactivator  recognition  surface  are 
inappropriately  positioned.  In  addition,  helix  12  itself  is  bound  such  that  it  partially  occludes  the 
static  region  of  the  coactivator-binding  surface  formed  by  residues  from  helices  3, 4,  and  5.  Yet, 
GEN  and  R,R-THC  clearly  exhibit  different  activities  on  ERP  in  mammalian  cells;  GEN  acts  as 
a  partial  agonist  (7)  and  R,R-THC  acts  as  a  pure  antagonist  (8, 9).  How  can  the  similarity  of  the 
two  structures  be  reconciled  with  the  different  activities  of  these  compounds  in  transcriptional 
assays?  The  simplest  model  that  explains  these  and  other  data  is  based  on  two  suppositions 
(Figure  5).  First,  helix  12  in  the  unliganded  ERP  LBD  is  in  equilibrium  between  the  R,R- 
THC/GEN-bound  conformation,  which  precludes  coactivator  association,  and  the  agonist-bound 
conformation  seen  in  the  E2-  and  DES-ERa  LBD  structures,  which  favors  coactivator 
association  (with  the  “inactive”  R,R-THC/GEN-bound  conformation  being  heavily  favored). 
Second,  rather  than  stabilizing  a  single  static  conformation  of  helix  12,  ligands  serve  to  shift  the 
balance  of  this  equilibrium.  It  is  also  highly  likely  that  the  coactivator  receptor-interaction 
domain  (e.g.  GRIP  NR  box  II  domain)  plays  an  active  role  in  shifting  this  equilibrium. 
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Apo 

inactive  active  +  coactivator 

inactive  ■■■■■ . .  active  full  agonist  (E2) 

inactive  «  r  active  partial  agonist  (GEN) 

inactive  active  antagonist  (THC) 

Prediction:  E2  >  GEN  >  Apo  >  THC 

Figure  5.  Equilibrium  model  of  Helix  12  positioning. 

R,R-THC  could  antagonize  receptor  activity  by  shifting  equilibrium  in  favor  of  the  inactive  R,R- 
THC/GEN-bound  conformation  (Figure  5).  If  R,R-THC  binding  forces  the  receptor  to  favor  this 
conformation  even  more  than  it  does  in  the  absence  of  ligand,  the  R,R-THC-receptor  complex 
would  be  expected  to  have  an  even  lower  affinity  for  coactivator  than  the  unliganded  receptor. 
This  would  be  consistent  with  the  absence  of  transcriptional  activity  of  the  R,R-THC-ER(3 
complex  and  the  conformation  of  helix  12  observed  in  the  crystal. 

Because  the  conformational  equilibrium  of  helix  12  would  be  difficult  to  observe  directly,  we 
chose  to  test  the  hypothesis  that  the  binding  of  different  ligands  should  modulate  the  affinity  of 
the  ERp  LBD  for  coactivator  such  that  E2-LBD  ~  DES-LBD  >  GEN-LBD  >  apo-LBD  >  R,R- 
THC-LBD.  The  affinities  of  various  ligand-ER|3  LBD  complexes  for  an  LXXLL  motif- 
containing  peptide  were  therefore  directly  measured  using  a  fluorescence  polarization-based 
binding  assay  (Appendix,  Figure  2b)  (3, 10).  The  full  agonist  complexes  bind  the  peptide 
approximately  1.5-fold  more  tightly  than  the  GEN  complex  (Kj[E2]  =  71  nM,  Kd[GEN]  =  104 
nM).  In  addition,  the  GEN  complex  binds  the  peptide  more  tightly  than  the  unliganded  receptor 
(Kd[apo]  =  215  nM),  whereas  the  R,R-THC  complex  (Ka[R,R-THC]  =  3.3  pM)  binds  the  peptide 
significantly  weaker  than  the  unliganded  receptor.  Thus,  these  data  are  entirely  consistent  with 
the  described  helix  12  conformational  equilibrium  model.  Furthermore,  recent  mutational  and 
crystallographic  data  suggest  that  the  positioning  of  helix  12  in  ERa  (11)  and  other  NRs  (12)  is 
also  dictated  by  a  ligand-sensitive  conformational  equilibrium. 

To  better  define  and  understand  “passive  antagonism”,  or  antagonism  without  side  chains,  one 
has  to  consider  the  differences  between  classical  estrogen  antagonists,  which  contain  bulky  side 
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chains,  and  molecules  like  R,R-THC.  Although  agonists  and  antagonists  bind  at  the  same  site 
within  the  core  of  the  LBD,  each  induces  distinct  conformations  in  the  transactivation  domain 
(AF-2)  of  the  LBD,  especially  in  the  positioning  of  helix  12,  providing  structural  evidence  for 
multiple  mechanisms  of  selective  antagonism  in  the  nuclear  receptor  family.  Interestingly,  the 
OHT/RAL  and  DES/E2  structures  collectively  reveal  and  define  a  multipurpose  docking  site  on 
ERa  that  can  accommodate  either  helix  12  or  one  of  several  coregulators.  In  addition,  a 
comparison  of  the  two  types  of  structures  reveals  that  there  are  at  least  two  distinct  mechanisms 
by  which  structural  features  of  OHT  promote  an  inhibitory  conformation  of  helix  12  (1, 6). 

Helix  12  positioning  is  determined  both  by  steric  considerations,  such  as  the  presence  of  an 
extended  side  chain  in  the  ligand,  and  by  local  structural  distortions  in  and  around  the  ligand 
binding  pocket.  Thus,  one  would  predict  that  effective  estrogen  antagonists  do  not  necessarily 
require  bulky  or  extended  side  chains. 

The  positioning  of  the  side  chains  of  OHT,  RAL,  and  ICI  directly  or  "actively"  preclude  the 
agonist-bound  conformation  of  helix  12  by  steric  hindrance.  Hence,  we  define  their  common 
mechanism  of  action  as  "active  antagonism".  Clearly,  RJR-THC  does  not  have  a  bulky  side 
chain  (Appendix,  Figure  la),  and  in  its  complex  with  ERf3,  helix  12  is  not  sterically  precluded 
from  adopting  the  agonist-bound  conformation  as  it  is  in  the  other  antagonist  complexes  (Figure 
2;  Appendix,  Figures  lb  and  2a).  Instead,  R,R-THC  antagonizes  ER(3  by  stabilizing  key  ligand 
binding  pocket  residues  in  nonproductive  conformations,  and  disfavoring  the  equilibrium  to  the 
agonist-bound  conformation  of  helix  12  (Figure  5).  Thus,  we  term  the  mechanism  of 
antagonism  of  R,R-THC  as  “passive  antagonism”.  Passive  antagonism  may  not  be  unique  to 
R,R-THC  and  ER|3.  There  are  many  examples  of  NR  ligands,  which  act  as  antagonists  even 
though  they  are  smaller  than  the  endogenous  agonists  of  these  NRs.  The  synthetic  androgen 
receptor  antagonist,  flutamide,  is  comparable  in  size  to  testosterone,  and  does  not  possess  an 
obvious  moiety  that  would  act  as  an  antagonist  side  chain  (13). 

Many  NRs  have  multiple  subtypes  that  possess  distinct  expression  patterns  and  that  regulate 
distinct  target  genes.  Antagonists  generated  through  the  addition  of  bulky  side  chains  to  agonist 
scaffolds  are  limited  to  being  antagonistic  on  one  or  more  subtypes  of  a  particular  NR.  The 
passive  antagonism  mechanism,  as  revealed  in  our  studies  through  direct  comparison  of  the  two 
R,R-THC-ER  LBD  complexes,  suggests  a  new  approach  to  achieving  NR  antagonism. 
Compounds  could  be  designed  to  selectively  stabilize  the  inactive  conformations  of  certain  NR 
subtypes  and  the  active  conformations  of  others.  Such  ligands  may  exert  novel  biological  and 
therapeutic  effects. 

A  manuscript  describing  the  above  data  was  recently  published  in  Nature  Structural  Biology  (3). 
Clearly,  significant  progress  was  made  during  the  past  three  years  in  our  efforts  to  understand  the 
unique  ER  subtype-selective  behavior  of  R,R-THC  and  the  R,R-THC-induced  conformation  of 
ER|3  that  contributes  to  “passive  antagonism”.  More  on  this  subject  is  included  in  the  summary 
of  Tasks  3  &  4  below. 

Tasks  3  &  4  (Third  year):  To  compare  and  correlate  structural  differences  among  hERa/{3- 
LBD  complexes  with  the  biological  activity  of  hERa/p. 
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We  have  made  significant  progress  in  trying  to  understand  the  contributions  of  residues  that  line 
the  ligand-binding  pocket  to  ligand  discrimination  between  ERa  and  ERfJ.  Progress  has  also 
been  made  in  our  effort  to  crystallize  ERa  bound  to  the  unique  SERM  GW5638.  However,  the 
main  effort  related  to  this  grant  during  the  third  year  was  to  better  understand  what  distinguishes 
ERa  from  ERp  in  their  ability  to  differentially  recognize  diverse  ligands  with  unique  behavior. 
To  achieve  this  goal,  we  generated  a  series  of  specific  mutant  ERs  and  tested  their  ability  to 
interact  with  known  coactivators  and/or  corepressors  and  to  act  as  selective  transcription 
modulators  in  co-transfection  assays.  Surprisingly,  we  discovered  that  residues  outside  the 
ligand-binding  pocket  play  important  roles  in  ligand  recognition  and  discrimination  between  the 
two  ER  subtypes. 

As  described  above,  R,R-THC  does  not  behave  as  an  ER|3  agonist.  Instead,  it  forms  alternative 
hydrogen  bonds  and  additional  packing  (van  der  Waals)  contacts  not  observed  in  the  ERa 
complex  (Appendix,  Figure  3b,  3c).  Hence,  R,R-THC  is  able  to  bind  to  the  ER|3  LBD  without 
nucleating  many  of  the  equivalent  cooperative  interactions  observed  in  the  ERa  complex.  For 
example,  in  the  ERa  complex  the  side  chain  of  Met  528  packs  against  the  side  chains  of  His  524 
and  Met  343,  stabilizing  the  helical  conformation  of  the  flanking  residues.  Because  R,R-THC 
binding  favors  alternative  conformations  of  His  475  and  Met  295  (equivalent  to  His  524  and  Met 
343,  respectively),  the  side  chain  of  Met  479  (equivalent  to  Met  528)  is  completely  disordered  in 
the  ERp  complex.  In  fact,  by  stabilizing  conformations  of  binding  pocket  residues  that  prevent 
their  interaction,  R,R-THC  binding  should  actually  disfavor  the  equilibrium  to  the  agonist-bound 
conformation  of  helix  12  (Figure  5;  Appendix,  Figure  2b). 


A  key  unanswered  question  relates  to  the  molecular/structural  differences  between  the  two  ER 
subtypes  that  allow  R,R-THC  to  stabilize  these  distinct  arrangements  of  ligand  binding  pocket 
residues.  There  are  only  two  binding  pocket  residues  which  are  not  identical  in  the  two 


R394/346 


Figure  6.  Superimposed  binding  pockets  of 
ERa  and  ER|3  bound  to  R,R-THC.  Brackets 
identify  residues  that  are  different  between  the 
two  ER  subtypes. 


subtypes:  Leu  384  (from  helix  5/6)  and  Met  421  (from  helix  8)  in  ERa,  which  are  equivalent  to 
Met  336  and  lie  373,  respectively,  in  ER(3  (Figure  6;  Appendix,  Figure  3c)  (14, 15).  These  two 
residue  pairs  are  similarly  positioned  on  opposite  faces  of  R,R-THC  in  their  respective 
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complexes.  In  the  ERa  complex,  Leu  384  and  Met  421  form  several  nonpolar  contacts  with 
R,R-THC  (Appendix,  Figure  3b).  By  contrast,  only  Met  336  interacts  with  R,R-THC  in  the 
ER|3  complex. 

The  two  ligand  binding  pocket  residues  that  differ  between  the  two  subtypes  may  partially 
explain  the  failure  of  R,R-THC  to  act  as  an  ERp  agonist.  Met  336  and  lie  373  in 
ERp  (equivalent  to  Leu  384  and  Met  421  in  ERa)  are  positioned  on  opposite  faces  of  R,R-THC 
(Appendix,  Figure  3b,  3c).  Modeling  (based  on  both  the  THC-  ERa  and  GEN-  ERp  structures) 
suggests  that,  because  of  the  different  sizes  of  the  side  chains  of  these  residues,  the  THC  A’  ring 
would  have  to  be  positioned  ~  2.4  A  lower  in  the  ERP  binding  pocket  relative  to  its  location  in 
the  crystal  in  order  for  THC  to  interact  with  Met  336  and  lie  373  and  to  orient  His  475  and  Leu 
476  in  productive  conformations  (data  not  shown).  This  lower  position  should  be  highly 
unfavorable  because  it  would  result  in  steric  clashes  between  the  B’  ring  ethyl  group  and  the  side 
chain  of  Leu  298  at  floor  of  the  binding  pocket  (Appendix,  Figure  3b,  3c).  R,R-THC 
presumably  binds  in  the  alternative  mode  observed  in  the  crystal  to  avoid  these  steric  clashes. 
Consistent  with  this  hypothesis,  the  dimethyl  analog  of  THC,  which  should  be  less  sterically 
hindered  from  adopting  the  lower  position,  acts  as  a  weak  ERP  partial  agonist  (8). 

Because  the  conformational  effects  of  ligand  binding  are  cooperative  in  nature,  residues  that  lie 
outside  the  binding  pocket  may  also  play  roles  the  subtype-selective  effects  of  R,R-THC. 
Nevertheless,  preliminary  mutagenesis  results  indicate  that  Leu  384/Met  336  and  Met  421/Ile 
373  do  contribute  to  these  effects  (Figure  7),  at  least  for  ERa.  However,  the  effect  is  only 
partial  for  ERa  and  the  corresponding  ERP  -*  ERa  amino  acid  swapping  mutations  (Met336  -» 
Leu384,  Ile373  ->  Met421)  do  not  result  in  the  recoveiy  of  any  detectable  R,R-THC  agonist 
activity  (Figure  7B).  The  same  result  was  observed  both  for  full  length  ERa  and  ERP  (data  not 
shown)  as  well  as  for  the  LBD  of  each  ER  fused  to  the  Gal4  DNA  binding  domain.  In  both 
cases,  a  transcriptional  activation  assay  was  used  to  assess  ER  or  ER  LBD  activity  in  transiently 
transfected  MDA-MB-231  breast  cancer  cells  in  which  recombinant  GRIP1  was  also  co¬ 
expressed  (see  attached  manuscript  for  experimental  details).  Estradiol  and  OHT  were  tested  as 
positive  and  negative  controls,  respectively  for  both  ER  subtypes  and  behaved  as  expected  with 
both  single  and  double  mutations  (Figure  7A).  In  total,  these  results  indicate  that  the  two  non- 
conserved  amino  acids  that  line  the  inside  of  the  LBD  pocket  (Figure  6)  are  not  sufficient  for 
selective  a  or  p  behavior  in  terms  of  response  to  R,R-THC.  However,  for  other  ligands,  such  as 
HPTE  (Figure  8),  these  two  amino  acids  do  appear  to  play  a  major  role  in  distinguishing 
between  ERa  and  ERp  behavior  (Figure  7B).  While  HPTE,  an  estrogenic  metabolite  of  the 
organochlorine  pesticide  methoxychlor,  is  a  potent  ERa  agonist  in  HepG2  cells  (16),  it  has 
minimal  agonist  activity  with  either  human  or  rat  ERp  and  almost  completely  abolishes 
estradiol-induced  ERp  activity.  Notably,  exchange  of  the  two  ERa/p-distinguishing  amino  acids 
effectively  converted  ERa  into  ERP  and  vice  versa  in  a  cell  based  reporter  assay  of  HPTE 
activity  (Figure  7B,  HPTE).  In  contrast,  the  behavior  of  certain  pyrazole  derivatives,  like 
propyl  pyrazole  triol  (Figures  7B  and  8,  PPT),  which  have  been  shown  to  be  selective  ERa 
agonists  with  little  or  no  affinity  for  ERp  (17),  are  unaffected  by  exchange  of  these  two  amino 
acids.  The  structural  basis  for  these  differences  has  yet  to  be  fully  determined,  although 
significant  progress  has  been  made  in  the  last  year,  as  described  below. 
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B  Gal4  ER  alpha  LBD  wt 

B  Gal4  ER  alpha  LBD 
L384M 
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M421I 

□  Gal4  ER  alpha  LBD 
L384M/M421 1 

B  Gal4  ER  beta  LBD  wt 


B  ER  alpha  LBD  wt 
BER  alpha  LBD  L384M 

□  ER  alpha  LBD  M421I 

□  ER  LBD  L384M/M421I 
b  ER  beta  LBD  wt 

B  ER  beta  LBD  M336L 
B  ER  beta  LBD  I373M 

□  ER  beta  LBD  M336L/I373M 


Figure  7.  Transcriptional  activity  of  Gal4-ERa/(3  LBD  mutants  versus 
wild  type  ERa/p  LBDs  (+  GRIP1;  in  MDA-MB23 1  cells).  A.  Positive 
and  negative  controls.  B.  ER  subtype-selective  ligands. 
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HPTE  (ERa  Agonists/ERp  Antagonists)  R,R-THC 

Figure  8.  Structures  of  ERa/p  subtype-selective  agonists  and  antagonists.  HPTE  = 
p-hydroxyphenyl-trichloroethane.  PPT  =  propyl  pyrazole  triol. 


Perhaps  these  contrasting  outcomes  are  not  quite  as  surprising  they  appear  since  we  and  others 
have  data  to  show  that  the  two  ER  subtypes  seem  to  prefer  different  default  conformations 
(Appendix,  Figure  2a)  and  that  even  for  the  same  ER  subtype,  structurally  similar  ligands  (e.g. 
OHT  and  RAL)  can  induce  unique  and  subtly  different  conformations.  Therefore,  it  is  very 
likely  that  amino  acids  outside  the  LBD  pocket  contribute  to  ligand  discrimination  and  ER 
subtype-selective  behavior.  Additional  proof  that  these  subtle  conformational  differences  are 
real  and  significant  comes  from  the  development  of  peptides  (generated  by  phage  display)  that 
can  distinguish  between  complexes  of  ERa  or  ERj3  bound  to  these  ligands  (18, 19).  For 
example,  a  peptide  designated  a/p  V  can  distinguish  between  the  OHT-ERa/p  complexes  and 
other  SERM-ER  complexes,  including  Ral-ERa/p  (18).  Other  peptides  (e.g.  7(5-16)  can 
selectively  recognize  GW5638-ERa/(5  compared  to  other  SERM-ER  complexes  (19).  One  of 
our  major  goals  is  to  understand  the  structural/molecular  basis  for  these  differences.  Thus,  one 
approach  is  to  determine  the  crystallographic  structures  of  ERa/ (3  LBD  complexes  with  these 
ligand-selective  peptides.  For  example,  we  are  currently  trying  to  crystallize  GW5638-ERa 
LBD  complexed  with  a  peptide  that  only  recognizes  this  ligand-ER  complex.  Interestingly,  this 
peptide,  which  was  identified  in  a  phage  display  screen,  has  a  high  degree  of  homology  with  the 
nuclear  receptor  interaction  domains  (CoRNR  boxes)  that  are  present  in  the  co-repressors  N-CoR 
and  SMRT  (20). 
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In  addition,  we  have  and  will  continue  to  exploit  site-directed  mutagenesis  of  both  ER  subtypes 
to  better  define  which  residues  distinguish  ERa  from  ERp  in  response  to  ligands  that  have 
unique  activity-  and/or  ER  subtype-selective  behavior.  An  example  of  the  value  of  this  approach 
is  shown  in  Figure  9.  By  introducing  additional  mutations  outside  the  residues  that  line  the 
ligand-binding  pocket  of  ER|3  LBD,  we  have  succeeded  in  recovering  partial  response  to  PPT, 
which  is  a  complete  ERa-selective  agonist.  The  choice  of  mutations  was  based  on  an  analysis  of 
differences  in  residue  interactions  between  the  two  ER  subtypes.  As  shown  in  Figure  10,  a 
series  of  hydrogen  bonds  (dashed  green  lines)  contribute  to  the  overall  structure  surrounding  the 
ligand-binding  pocket  of  ERa.  These  interactions  are  not  present  in  ER(3.  We  therefore  mutated 
the  corresponding  ER(3  residues  to  re-create  these  essential  interactions.  Interestingly,  a  cluster 
of  five  mutations  (Mut  5)  is  more  effective  than  adding  two  more  mutations  to  the  same  cluster 


Figure  9.  Transcriptional  activity  of  ER(3  mutants  in  response  to  PPT  (ERE- 
tk-Luc;  Ishikawa  Cells).  Dose-response  curves.  Beta  =  wild  type  ERfJ;  Mut 
2  =  M336L+ I373M;  Mut  5  =  Mut  2  +  M296G  +  L300N  +  L304R;  Mut  7  = 
Mut  5  +  I282L  +  D359N.  E2  =  estradiol-mediated  activity.  Y  axis  =  fold 
induction.  X  axis  =  log  dose  [ligand].  VEH  =  vehicle. 


(Mut  7)  (Figure  9),  demonstrating  that  not  all  ERa/p  residue  exchanges  are  additive  in 
effectiveness.  The  divergence  of  the  two  ER  subtypes  undoubtedly  took  a  long  time  to  evolve, 
which  is  why  we  are  trying  to  use  a  rational  structure-based  approach  coupled  with  functional 
assays  to  understand  the  nature  of  these  differences.  Despite  considerable  effort,  we  have  yet  to 
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determine  the  minimum  number  of  residues  that  are  required  for  ER|3  to  recognize  R,R-THC  as 
an  agonist.  Also,  the  partial  recovery  of  PPT  agonism  observed  in  the  Mut  5  ER(3  (Figure  9) 
demonstrates  that  pocket  and  surface  residue  requirements  are  not  the  same  for  different  ligands 
and  the  relationship  between  the  LBDs  of  the  two  ER  subtypes  is  more  complicated  than  we  first 
imagined.  These  data  clearly  illustrate  that  secondary  structure  interactions  and  flexibility  are 
key  determinants  of  ligand  specificity  for  the  two  ER  subtypes,  revealing  novel  aspects  of 
nuclear  receptor  function.  Equally  important,  residues  outside  the  ligand-binding  pocket  appear 
to  be  critical  for  ligand  discrimination  between  two  closely  related  nuclear  receptors. 


Figure  10.  Hydrogen  bonds  that  are  ERa-specific.  Helix  12  is  shown  at  the  far 
left  in  this  partial  view  of  the  R,R-THC-ERa  LBD  structure.  The  GRIP  peptide 
(upper  left  comer)  is  colored  red. 


We  will  continue  to  use  a  structure-based  mutagenic  approach  to  identify  key  residues  that  are 
responsible  for  ligand  discrimination  between  the  two  ER  subtypes.  This  type  of  information 
will  be  essential  for  our  compound  modeling  and  drug  design  efforts.  We  will  want  to  know,  for 
example,  the  minimum  number  of  amino  acids  required  for  ER  subtype  selection  for  each  of  the 
compounds  mentioned  above.  We  will  also  need  to  know  more  about  the  conformational 
flexibility  of  the  two  ER  subtypes  and  how  that  flexibility  distinguishes  between  coactivator  and 
corepressor  preferences.  The  success  of  this  project  has  therefore  not  only  been  very  informative 
in  its  own  right  but  has  also  led  to  new  approaches  that  will  help  in  the  design  of  molecules  that 
may  have  ER  subtype-selective  therapeutic  applications  in  HRT  and  breast  cancer  management 
or  prevention. 
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KEY  RESEARCH  ACCOMPLISHMENTS 

•  Determination  and  refinement  of  crystallographic  structures  of  ERa  and  ER|3  LBDs  bound  to 
R,R-THC,  an  ERa  agonist/ER(3  antagonist.  New  insights  into  mechanisms  of  agonism  and 
antagonism  revealed. 

•  Molecular  basis  of  novel  concept  “passive  antagonism”,  or  antagonism  without  side  chains, 
predicted  from  structures  of  ERa  and  ER(3  LBDs  bound  to  R,R-THC  and  TIF2/GRIP1  NR 
box  II  peptide. 

•  Novel  concept  of  helix  12  dynamic  equilibrium:  Positioning  of  helix  12  is  more  variable  than 
previously  believed  and  is  a  key  discriminator  among  diverse  agonist-  and  antagonist- 
induced  conformations  of  ERa  and  ER(3.  Hydrophobic  cleft  in  the  LBD  can  accommodate 
coactivators  (agonists),  corepressors  (antagonists),  and/or  different  orientations  of  helix  12  in 
the  presence  of  different  antagonists. 

•  Transcriptionally  active  conformation  of  ER  LBD  requires  cooperative  participation  of 
ligand,  LBD,  and  coactivator. 

•  Molecular  basis  of  ER  subtype-selective  agonism  and  antagonism  predicted  from  structure  of 
ERa  and  ER(3  LBDs  bound  to  R,R-THC.  R,R-  THC  binds  to  the  ER(3  LBD  without 
nucleating  many  of  the  equivalent  cooperative  interactions  observed  in  the  ERa  complex  that 
reflect  agonist  behavior. 

•  Mutagenic  analysis  of  ERa  and  ER|3  LBDs  reveals  that  ligand-selective  a  and  (5  character, 
especially  transcriptional  activity,  depend  on  more  than  the  two  amino  acids  that  differ  in  the 
ligand  binding  pockets  of  the  two  ERs. 

•  Amino  acids  that  are  required  for  PPT  (a  potent  ERa-selective  agonist)  recognition  by  ER|3 
are  identified  by  mutagenic  analysis  combined  with  functional  assays.  These  resides  are  not 
sufficient  for  ERp  recognition  of  R,R-THC  as  an  agonist. 

•  LBD  pocket  and  surface  residue  requirements  are  not  the  same  for  different  ligands  and  the 
relationship  between  the  LBDs  of  the  two  ER  subtypes  is  more  complicated  than  first 
imagined.  Secondary  structure  interactions  and  flexibility  are  key  determinants  of  ligand 
specificity  for  the  two  ER  subtypes,  revealing  novel  aspects  of  nuclear  receptor  function 

REPORTABLE  OUTCOMES 

•  Manuscripts,  Abstracts,  Presentations: 

A  manuscript  detailing  the  results  reported  here  was  recently  published  in  Nature 

Structural  Biology  (3).  It  is  included  in  the  Appendix. 

Abstracts  for  two  meetings  are  included  in  the  Appendix. 

Presentations: 

Seminar  Speaker,  Molecular  and  Cellular  Biology  Research  Seminar  Series,  Sunnybrook 
and  Women’s  College  Health  Sciences  Centre,  Toronto,  Ontario,  Canada,  11/01 

Symposium  Speaker,  Workshop  on  Genomic  vs  Non-Genomic  Steroid  Actions: 

Encountered  or  Unified  Views,  Juan  March  Institute,  Madrid,  Spain,  12/17-19, 2001. 

Symposium  Speaker,  Breast  Cancer  Symposium:  Think  Tank  12,  St.  Maarten,  The 
Netherlands  Antilles,  1/02 

Symposium  Speaker,  Sixth  Annual  Breast  Cancer  Symposium,  Sun  Valley,  Idaho,  3/02 

Speaker,  Sixth  Annual  Workshop,  Frontiers  in  Estrogen  Action,  Amelia  Island,  FL,  4/02 

•  Patents:  none 

•  Degrees  obtained  with  award  support:  none 


17 


Geoffrey  L.  Greene,  Ph.D. 


•  Development  of  cell  lines,  etc:  none 

•  Informatics: 

R,R-THC-ERa/p  crystallographic  coordinates  have  been  deposited  in  the  PDB  and  were 
released  upon  publication  of  the  structures. 

•  Funding  applied  for  based  on  work  supported  by  this  award: 

An  R01  proposal  (CA89489),  "Development  and  Characterization  of  Novel  SERMs” 
was  funded  in  August,  2001.  This  project  will  continue  and  extend  the  structural 
work  supported  by  this  DOD  award.  In  addition,  a  drug  discovery  project, 

“Development  of  Novel  Reagents  for  Breast  Cancer  Prevention  and  Treatment”, 
which  complements  this  DOD  award  was  funded  January  1, 2002  by  the  Ludwig 
Fund  for  Cancer  Research. 

•  Employment  or  research  opportunities  applied  for:  none 
CONCLUSIONS 

The  major  goal  of  this  investigation  was  to  determine  the  three-dimensional  structures  of  the  two 
known  human  estrogen  receptors  (ERa  and  ER|3)  complexed  with  receptor-selective  estrogens 
and  antiestrogens  (SERMs).  Ligands  that  act  differentially  on  ERa  and  ERf3  are  potentially 
valuable  both  as  tools  to  dissect  the  biological  roles  of  the  two  ERs  and  as  novel  therapeutics 
with  pharmacological  properties  distinct  from  existing  drugs.  The  synthetic  compound,  R,R- 
5,1  l-cw-diethyl-5, 6,11, 12-tetrahydrochry sene-2, 8-diol  (R,R-THC),  was  previously  identified  as 
a  selective  estrogen  agonist  when  bound  to  ERa  and  as  an  antagonist  when  bound  to  ER(3.  To 
better  understand  this  selective  behavior,  the  crystallographic  structures  of  human  ERa  and  ERfl 
ligand-binding  domains  (LBDs)  complexed  with  R,R-THC  were  solved.  Consistent  with  the 
prediction  that  bulky/extended  side  chains  are  not  essential  for  antagonist  behavior,  R,R-THC 
antagonizes  ER(J  in  a  manner  very  different  from  4-hydroxytamoxifene  (OHT),  raloxifene 
(RAL)  and  ICI  182,780  (ICI).  The  positioning  of  the  side  chains  of  OHT,  RAL,  and  ICI  directly 
or  "actively"  preclude  the  agonist-bound  conformation  of  helix  12  by  steric  hindrance.  Hence, 
we  define  their  common  mechanism  of  action  as  "active  antagonism".  R,R-THC  does  not  have 
an  extended  side  chain  and  in  its  complex  with  ER|3,  helix  12  is  not  sterically  precluded  from 
adopting  the  agonist-bound  conformation  as  it  is  in  the  other  antagonist  complexes.  Instead, 
R,R-THC  antagonizes  ERp  by  stabilizing  key  ligand-binding  pocket  residues  in  noninteracting 
conformations,  and  disfavoring  the  equilibrium  to  the  agonist-bound  conformation  of  helix  12. 
Thus,  we  term  the  mechanism  of  antagonism  by  R,R-THC  as  “passive  antagonism”.  This 
mechanism  suggests  a  new  approach  to  achieving  NR  antagonism.  Compounds  could  be 
designed  to  selectively  stabilize  the  inactive  conformations  of  certain  NR  subtypes  and  the  active 
conformations  of  others,  thereby  eliciting  novel  biological  and  therapeutic  effects. 

In  order  to  determine  which  amino  acids  are  responsible  for  ligand  discrimination  between  ERa 
and  ERp,  we  mutagenized  the  LBDs  of  both  ERs,  starting  with  the  two  residues  (L384/M336, 
M421/I373)  that  differ  in  the  pockets  of  each  ER  subtype.  PPT  (propyl  p>yrazole  triol),  an  ERa- 
selective  agonist  with  very  low  affinity  for  ER(5,  as  well  as  HPTE  (2,2-bis(p-hydroxyphenyl)-l,l, 
1-trichloroethane)  and  THC,  both  of  which  act  as  ERa  agonists  and  ER(3  antagonists,  were  used 
to  assess  the  consequences  of  these  mutations.  Interestingly,  conversion  of  the  two  pocket 
residues  in  ERa  to  the  corresponding  residues  in  ERj3  markedly  reduced  the  activity  of  HFTE 
and  THC  on  ERa,  but  had  no  effect  on  PPT  activity.  The  converse  mutations  in  ERp  afforded 
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moderate  agonist  activity  with  HPTE  but  no  activity  with  either  PPT  or  THC.  Thus,  pocket  and 
surface  residue  requirements  are  not  the  same  for  different  ligands  and  the  relationship  between 
the  LBDs  of  the  two  ER  subtypes  is  more  complicated  than  first  imagined.  The  crystal  structures 
of  THC-ERa  and  ER(5  LBDs  suggested  additional  residues  that  might  contribute  to  ligand 
selectivity  through  altered  secondary  structure  interactions.  Mutation  of  some  of  these  residues 
in  ERp  resulted  in  partial  recovery  of  PPT  agonist  activity  but  were  not  sufficient  to  convert 
R,R-THC  from  an  ER0  antagonist  to  an  agonist.  These  data  illustrate  that  secondary  structure 
interactions  and  flexibility  are  key  determinants  of  ligand  specificity  for  the  two  ER  subtypes, 
revealing  novel  aspects  of  nuclear  receptor  function.  Equally  important,  residues  outside  the 
ligand-binding  pocket  appear  to  be  critical  for  ligand  discrimination  between  two  closely  related 
nuclear  receptors. 

So  what?  It  is  clear  from  the  various  solved  structures  for  ERa  and  ER|3  bound  to  diverse 
SERMs  that  the  resulting  ER  conformations  and  ligand  pharmacologies  are  not  as  simple  as 
previously  believed.  Rather  than  just  two  distinct  conformations  of  ER  that  reflect  agonism  or 
antagonism,  it  appears  that  multiple  conformations  are  possible,  each  of  which  may  mediate  a 
different  pharmacology  for  a  particular  ligand  in  a  given  tissue,  such  as  the  breast,  uterus,  bone, 
or  cardiovascular  system.  It  has  become  clear  that  variable  antagonism  can  be  achieved  both  by 
the  presence  of  bulky  side  chains  and  by  molecules  that  lack  bulky  side  chains.  Thus,  it  is 
important  to  solve  additional  structures  for  ERa  and  ERp  complexed  with  ligands  that  have  been 
shown  to  have  different  tissue-selective  pharmacologies,  such  as  GW5638  and  ICI  182,780 
(Faslodex).  In  addition,  structure-based  mutagenic  analyses  of  the  two  ER  subtypes  will  help 
define  residues  that  contribute  to  a  and  (3  activity-selective  and  subtype-selective  behavior.  Such 
information  is  essential  to  design  compounds  that  have  defined  properties.  A  major  goal  of  this 
investigation  is  to  characterize  and/or  design  compounds  that  are  useful  as  chemopreventive 
agents,  especially  for  breast  and  uterine  cancers.  To  be  useful  in  a  clinical  setting,  such  as 
hormone  replacement  therapy,  such  compounds  must  also  have  the  beneficial  effects  of  estrogen, 
including  maintenance  of  bone  density,  protection  of  the  cardiovasculature,  prevention  of  hot 
flashes  and  perhaps  delayed  onset  of  Alzheimer’s  disease  and  maintenance  of  cognitive  function. 
It  seems  unlikely  that  a  single  molecule  will  have  all  of  these  properties.  However,  it  is  clear  that 
improved  SERMs  can  be  designed  and/or  discovered.  Since  most,  if  not  all,  of  the  behavior  of 
the  two  known  ERs  is  dictated  by  the  conformations  induced  by  diverse  natural  and  synthetic 
ligands,  determination  and  analysis  of  the  corresponding  3D  structures  for  representative  SERMs 
should  be  a  valuable  tool  for  dissecting  underlying  molecular  mechanisms  and  for  designing  new 
generation  compounds.  We  will  therefore  continue  our  efforts  to  determine  and  correlate  ER- 
SERM  structures  with  their  known  biological  activities  and  to  understand  the  unique 
structural/molecular  characteristics  that  define  and  differentiate  ERa  and  ERj3.  These  approaches 
should  help  in  the  design  of  molecules  that  may  have  ER  subtype-selective  therapeutic  or 
preventive  applications  in  HRT,  hormone/tissue-specific  diseases  and  breast  cancer. 
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Workshop  on  Genomic  vs  Non-Genomic  Steroid  Actions:  Encountered  or  Unified  Views,  Juan 
March  Institute,  Madrid,  Spain,  12/17-19,  2001. 

ERa  and  ER(3:  Complexities  of  ligand  recognition  and  its  consequences 

Geoffrey  L.  Greene.  The  Ben  May  Institute  for  Cancer  Research,  The  University  of  Chicago, 
Chicago,  Illinois  60637 

Estrogens  and  SERMs  (selective  estrogen  receptor  modulators)  regulate  diverse  cellular 
activities  via  one  or  both  of  two  known  estrogen  receptor  subtypes  (ERa  and  ER(3)  in  hormone 
responsive  tissues  and  cancers.  Liganded  ERs  can  interact  with  a  complex  mix  of  coactivators, 
corepressors  and  other  signaling  molecules  that  differ  in  expression  and  importance  from  tissue 
to  tissue.  In  addition,  different  SERMs  may  alter  the  affinity  and/or  selectivity  of  one  or  both 
ERs  for  these  coregulators,  allowing  for  tissue  selective  responses.  Recently,  insight  into  the 
molecular  basis  of  estrogen  agonism  and  antagonism  has  been  revealed  by  the  crystal  structures 
of  ERa  and  ER|3  ligand  binding  domains  (LBDs)  complexed  with  several  ligands,  including 
estradiol  (E2),  diethylstilbestrol  (DES),  raloxifene  (RAL),  4-hydroxytamoxifen  (OHT),  and  the 
phytoestrogen  genistein  (GEN).  For  agonists  like  DES,  inclusion  of  a  peptide  derived  from  an 
essential  LXXLL  interaction  motif  (NR  box)  found  in  several  related  pl60  nuclear  receptor 
transcriptional  co-activators  has  helped  define  the  AF-2/co-activator  interface  on  ER.  Although 
agonists  and  antagonists  bind  at  the  same  site  within  the  core  of  the  LBD,  each  induces  distinct 
conformations  in  the  transactivation  domain  (AF-2)  of  the  LBD,  especially  in  the  positioning  of 
helix  12,  providing  structural  evidence  for  multiple  mechanisms  of  selective  antagonism  in  the 
nuclear  receptor  family.  Interestingly,  the  OHT/RAL  and  DES/E2  structures  collectively  reveal 
and  define  a  multipurpose  docking  site  on  ERa  that  can  accommodate  either  helix  12  or  one  of 
several  coregulators.  In  addition,  a  comparison  of  the  two  structures  reveals  that  there  are  at 
least  two  distinct  mechanisms  by  which  structural  features  of  OHT  promote  an  inhibitory 
conformation  of  helix  12.  Helix  12  positioning  is  determined  both  by  steric  considerations,  such 
as  the  presence  of  an  extended  side  chain  in  the  ligand,  and  by  local  structural  distortions  in  and 
around  the  ligand  binding  pocket.  Thus,  one  would  predict  that  effective  estrogen  antagonists  do 
not  necessarily  require  bulky  or  extended  side  chains. 

Ligands  that  act  differentially  on  ERa  and  ER(3  are  potentially  valuable  both  as  tools  to  dissect 
the  biological  roles  of  the  two  ERs  and  as  novel  therapeutics  with  pharmacological  properties 
distinct  from  existing  drugs.  As  part  of  a  search  for  ER  subtype-selective  ligands,  the  synthetic 
compound,  R,R-5,ll-ds-diethyl-5,6,ll,12-tetrahydrochrysene-2,8-diol  (R,R-THC),  was 
identified  as  a  selective  estrogen  agonist  when  bound  to  ERa  and  as  an  antagonist  when  bound 
to  ER|3.  To  better  understand  this  selective  behavior,  the  crystallographic  structures  of  human 
ERa  and  ER(3  ligand  binding  domains  (LBDs)  complexed  with  R,R-THC  were  solved  and 
refined,  suggesting  mechanisms  by  which  this  compound  can  act  as  an  ERa  agonist  and  as  an 
ER|3  antagonist.  Consistent  with  the  prediction  that  bulky/extended  side  chains  are  not  essential 
for  antagonist  behavior,  R,R-THC  antagonizes  ER(3  in  a  manner  very  different  from  OHT  and 
RAL.  The  positioning  of  the  side  chains  of  OHT,  RAL,  and  ICI  directly  or  "actively"  preclude 
the  agonist-bound  conformation  of  helix  12  by  steric  hindrance.  Hence,  we  define  their  common 
mechanism  of  action  as  "active  antagonism".  R,R-THC  does  not  have  a  bulky  side  chain  and  in 
its  complex  with  ER(3,  helix  12  is  not  sterically  precluded  from  adopting  the  agonist-bound 
conformation  as  it  is  in  the  other  antagonist  complexes.  Instead,  R,R-THC  antagonizes  ER|3  by 
stabilizing  key  ligand  binding  pocket  residues  in  noninteracting  conformations,  and  disfavoring 
the  equilibrium  to  the  agonist-bound  conformation  of  helix  12.  Thus,  we  term  the  mechanism  of 
antagonism  of  R,R-THC  as  “passive  antagonism”.  The  passive  antagonism  mechanism,  as 
revealed  here  through  direct  comparison  of  the  two  R,R-THC-ER  LBD  complexes,  suggests  a 
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Molecular  Determinants  of  Ligand  Selectivity  Between  ERa  and  ERp 
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The  development  of  ligands  that  are  selective  for  ERa  or  ERp  provides  an  important  tool 
for  dissecting  their  different  roles  in  vivo,  and  for  understanding  the  structural  basis  of 
their  ER  subtype  selectivity.  In  order  to  delineate  which  amino  acids  determine  ligand 
selectivity,  we  undertook  a  mutational  analysis  within  the  ligand-binding  domains  of  the 
full-length  estrogen  receptors,  starting  with  the  two  residues  (L384/M336,  M42 1/1373) 
that  differ  in  the  ligand-binding  pockets  of  each  ER  subtype.  We  examined  PPT  (propyl 
pyrazole  triol),  an  ERa-selective  agonist  with  very  low  affinity  for  ERp,  as  well  as  HPTE 
(2, 2-bis(p-hydroxy phenyl)- 1,1, 1-trichloroethane)  andTHC  (R,R-5,1 1-c/s-diethyl- 
5,6,1 1,12-  tetrahydro-chrysene-2,8-diol),  which  act  as  agonists  on  ERa  but  as  antagonists 
on  ERp.  Mutation  of  the  two  ligand-binding  pocket  residues  in  ERa  to  the  corresponding 
residues  in  ERp  reduced  the  activity  of  HPTE  and  THC  on  ERa,  but  had  no  effect  on 
PPT  activity.  The  converse  mutations  in  ERp  produced  moderate  agonist  activity  only 
with  HPTE.  The  recently  solved  crystal  structures  of  THC  bound  to  the  ERa  and  ERp 
LBDs  suggested  additional  residues  on  the  surface  that  might  contribute  to  ligand 
selectivity  through  altered  secondary  structure  interactions.  THC  does  not  enter  into  the 
pocket  as  deeply  in  ERp  as  in  ERa,  which  results  in  a  flexing  of  helix  1 1  away  from  the 
pocket  to  accommodate  the  ligand.  Both  the  flexibility  of  helix  1 1  and  the  position  of 
helix  12  can  be  explained  by  amino  acid  differences  in  helix  3.  ERa  Gly  344  and  Asn 
348  stabilize  the  agonist  position  by  interacting  with  the  loop  between  helix  1 1  and  12, 
and  forming  a  hydrogen  bond  with  the  beginning  of  helix  12,  respectively.  The 
corresponding  ERp  Met  296  and  Lys  300  residues  do  not  promote  these  agonist 
stabilizing  interactions,  suggesting  that  helix  12  may  be  preferentially  positioned  away 
from  helix  1 1  to  allow  greater  flexibility  in  accommodating  ligands.  This  flexibility  is 
needed  because  the  back  of  the  ligand  pocket  appears  to  be  narrower  in  ERp.  The 
distance  between  the  beta  sheet  and  helix  3  is  linked  with  six  additional  hydrogen  bonds 
in  ERa  that  are  not  present  in  ERp,  associated  with  residues  Ser329,  Asn352,  and 
Asn407.  Mutational  analysis  revealed  that  PPT  and  HPTE  act  as  partial  agonists  on  ERp 
after  mutating  the  ligand  pocket  and  surface  residues  described  above  to  the 
corresponding  residues  in  ERa.  HPTE  displayed  less  agonist  activity  with  the  ligand 
pocket  mutants  alone.  In  contrast,  additional  mutations  appear  to  be  required  for  ERp  to 
recognize  PPT  as  an  agonist,  suggesting  that  ligands  differ  in  the  relative  contributions 
that  ligand  pocket  and  surface  residues  make  in  determining  receptor  activity.  This  work 
illustrates  that  secondary  structure  interactions  and  flexibility  are  key  determinants  of 
ligand  specificity  for  the  two  ER  subtypes,  revealing  novel  aspects  of  nuclear  receptor 
function.  Equally  important,  this  may  be  the  first  example  of  nuclear  receptor  subtype 
specificity  with  determinants  outside  the  ligand-binding  pocket. 
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The  R,R  enantiomer  of  5,ll-cis-diethyl-5,6,ll,12-tetra- 
hydrochry$ene-2,8-diol  (THC)  exerts  opposite  effects  on  the 
transcriptional  activity  of  the  two  estrogen  receptor  (ER) 
subtypes,  ERa  and  ERp.  THC  acts  as  an  ERa  agonist  and  as 
an  ERp  antagonist.  We  have  determined  the  crystal  struc¬ 
tures  of  the  ERa  ligand  binding  domain  (LBD)  bound  to 
both  THC  and  a  fragment  of  the  transcriptional  coactivator 
GRIP1,  and  the  ERp  LBD  bound  to  THC.  THC  stabilizes  a 
conformation  of  the  ERa  LBD  that  permits  coactivator 
association  and  a  conformation  of  the  ERp  LBD  that 
prevents  coactivator  association.  A  comparison  of  the  two 
structures,  taken  together  with  functional  data,  reveals  that 
THC  does  not  act  on  ERp  through  the  same  mechanisms 
used  by  other  known  ER  antagonists.  Instead,  THC  antago¬ 
nizes  ERp  through  a  novel  mechanism  we  term  passive 
antagonism’. 

The  physiological  effects  of  both  endogenous  and  synthetic 
estrogens  are  mediated  by  the  estrogen  receptors  (ERs),  ERa  and 
ERp,  which  are  members  of  the  nuclear  receptor  (NR)  super¬ 
family  of  ligand-regulated  transcription  factors1,2.  In  a  search  for 
ER  subtype-selective  ligands,  the  R,R  enantiomer  of  5,11-ds- 
diethyl-5,6,ll,12-tetrahydrochrysene-2,8-diol  (THC)  (Fig.  la) 
was  determined  to  be  a  novel  ER  ligand,  with  an  approximate 
six-fold  affinity  preference  for  ERp  over  ERa  and  profoundly 
different  activation  profiles  on  the  two  ERs3,4.  THC  functions  as 
an  agonist  on  ERa,  with  an  EC50  of  ~3  nM,  but  has  no  effect  on 
the  transcriptional  activity  of  ERp.  Instead,  this  compound 
potently  antagonizes  the  effects  of  the  endogenous  estrogen 
17p-estradiol  (E2)  on  ERp,  with  an  IC50  <  10  nM.  Consistent 
with  these  subtype-specific  effects  on  transcriptional  activation, 
THC  also  affects  the  two  ERs  differently  in  protease  protection 
and  other  biochemical  assays5. 

The  C-terminal  ligand-binding  domains  (LBDs)  of  ERa  and 
ERP  each  contain  a  ligand-responsive  transcriptional  activation 
function  (AF-2)6-9.  Agonists,  such  as  E2  and  diethylstilbestrol 


Fig.  1  Overall  structures  of  the  THC-ER  LBD  complexes,  a,  Chemical 
structures  of  17p-estradioi  (E2),  ICI  164,384  (ICI),  diethylstilbestrol  (DES), 
4-hydroxytamoxifen  (OHT)  and  /?,rt-5,11-e/s-diethyl-5,6,1 1,12-tetrahy- 
drochrysene-2,8-dio!  (THC).  The  side  chains  of  ICI  and  OHT  are  high¬ 
lighted  in  orange,  b,  Two  equivalent  orthogonal  views  of  the  THC-ERa 
LBD-GR1P1  NR  box  II  peptide  and  the  THC-ERp  LBD  complexes,  showing 
that  the  two  ER  LBDs  adopt  distinct  conformations  when  bound  to  THC. 
Both  the  ERa  and  ER(3  LBDs  are  depicted  in  ribbon  representation  and 
colored  light  green  and  blue,  respectively.  In  both  complexes,  helix  12  is 
colored  red,  and  THC  (green)  is  shown  in  space-filling  representation.  In 
the  ERa  complex,  the  coactivator  peptide  is  depicted  as  a  purple  ribbon, 
and  helices  3,  4,  5  and  12  have  been  labeled  H3,  H4,  H5,  and  H12.  Dashed 
lines  represent  regions  of  the  structures  that  have  not  been  modeled. 
Panel  ( b )  and  Figs  2a,  3a, c  and  4  were  generated  using  BOBSCRIPT32  and 
rendered  using  Raster3D33. 


(DES),  increase  AF-2  activity.  In  contrast,  pure  antagonists,  such 
as  ICI  164,384  (ICI),  and  mixed-function  partial  agonist/antag¬ 
onists,  such  as  4-hydroxytamoxifen  (OHT)  and  raloxifene 
(RAL),  block  AF-2  activity10-12  (Fig.  la).  The  effects  of  THC  on 
the  transcriptional  activity  of  fusions  between  the  GAL4  DNA- 
binding  domain  and  the  two  ER  LBDs  indicate  that  THC  stimu¬ 
lates  ERa  AF-2  activity  and  inhibits  agonist-induced  ER|3  AF-2 
activity  (data  not  shown),  consistent  with  its  effects  on  the  full- 
length  receptors. 

Structural  studies  indicate  that  NR  ligands  regulate  AF-2 
activity  by  modulating  the  structure  of  NR  LBDs13’14.  NR  ligands 
bind  to  their  cognate  receptors  within  a  hydrophobic  pocket 
formed  within  the  core  of  the  narrower  or  'lower5  half  of  the 
LBD.  The  binding  of  an  agonist  to  an  NR  LBD  stabilizes  the  posi¬ 
tioning  of  the  most  C-terminal  helix  of  the  LBD,  helix  12,  such 
that  it  lies  across  the  opening  to  the  binding  pocket  formed  by 
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Fig.  2  Conformational  equilibrium  of  helix  12. 
a,  Superpositions  of  the  Ca  trace  of  the 
THC-ERp  LBD  complex  (light  blue)  on  those  of 
the  GEN-ERp  LBD  (yellow)  and  the  RAL-ERp 
LBD  (orange)  complexes.  Helix  12  in  the  RAL 
and  GEN  structures  is  colored  red.  Least  squares 
superpositions  were  generated  using  LSQ- 
MAN34  (THC/RAL  has  1.0  A  r.m.s.  deviation  over 
196  matched  atoms  with  a  3.8  A  cutoff,  and 
THC/GEN  has  1.1  A  r.m.s.  deviation  over  218 
matched  atoms  with  a  3.8  A  cutoff).  Helix  12 
from  each  of  the  monomers  in  the  THC-ERp, 
GEN-ERp  and  THC-ERa  crystals  is  located  in  a 
distinct  packing  environment,  indicating  that 
the  similarities/differences  in  helix  12  position¬ 
ing  among  these  complexes  are  not  the  conse¬ 
quence  of  crystal  lattice  effects,  b,  Equilibrium 
binding  of  a  rhodamine-labeled  Leu-X-X-Leu- 
Leu  motif-containing  peptide  to  the  ERp  LBD 
alone  or  bound  to  E2,  DES,  GEN  and  THC  was 
analyzed  by  measuring  fluorescence  polariza¬ 
tion  as  a  function  of  receptor  or  ligand-recep¬ 
tor  complex  concentration.  Values  represent 
the  mean  ±  s.d.  for  each  experimental  condi¬ 
tion  performed  in  triplicate,  and  the  data  were 
fit  using  nonlinear  least  squares  analysis  (Kd[E2] 
=  71  ±  2  nM,  Kd[GEN]  =  104  ±  4  nM,  Kd[apo]  = 
215  ±  5  nM  and  Kd[THC]  =  3.3  ±  0.3  pM). 


helices  3,  5/6  and  11.  This  conformation  of  the  LBD  allows  the 
receptor  to  interact  with  a  class  of  proteins  known  as  transcrip¬ 
tional  coactivators,  which  mediate  ligand-dependent  transcrip¬ 
tion  of  NRs15.  pi 60  coactivators  recognize  agonist-bound  NR 
LBDs  via  the  short  sequence  motif,  Leu-X-X-Leu-Leu  (where  X 
is  any  amino  acid),  called  the  NR  box15.  NR  boxes  form  amphi- 
pathic  a-helices  that  recognize  a  hydrophobic  groove  on  the  sur¬ 
face  of  an  agonist-bound  LBD  formed  by  residues  from  helices  3, 
4,  5  and  12  (ref.  13).  NR  AF-2  antagonists  sterically  preclude  the 
proper  formation  of  the  NR  box-binding  site13,14.  The  structures 
of  the  BMS614-retinoic  acid  receptor~aK\  OHT-ERa17, 
RAL-ERa18,  RAL-ERp19  and  ICI-ERP20  LBD  complexes  reveal 
that  each  of  these  antagonists  possesses  a  bulky  side  chain 
(Fig.  la)  that  cannot  be  contained  within  the  ligand  binding 
pocket.  These  side  chains  protrude  out  of  the  opening  to  the 
binding  pocket  formed  by  helices  3, 5/6  and  1 1 ,  thereby  prevent¬ 
ing  helix  12  from  adopting  the  agonist-bound  conformation. 
The  positioning  of  the  side  chains  of  BMS614,  OHT  and  RAL 
forces  helix  12  to  adopt  an  alternative  conformation  in  which  it 
binds  to  and  occludes  the  remainder  of  the  coactivator  binding 
site  by  mimicking  the  interactions  formed  by  the  NR  box16-19. 
The  ICI  side  chain,  in  contrast,  binds  directly  to  the  coactivator¬ 
binding  site  of  ERp,  causing  helix  12  to  be  completely  disor¬ 
dered20.  The  structure  of  the  ERP  LBD  bound  to  the  ERp  partial 
agonist  genistein  (GEN)  reveals  that  ligand  binding  can  stabilize 
yet  another  conformation  of  helix  12  (ref.  19).  In  this  complex, 
helix  12  is  bound  over  the  ligand-binding  pocket  in  a  position 
such  that  it  occludes  the  coactivator  recognition  surface  only 
partially.  However,  the  functional  significance  of  this  conforma¬ 
tion  of  the  LBD  is  unclear. 

Given  its  behavior  in  transactivation  and  other  assays3-5, 
THC  would  be  predicted  to  both  stabilize  a  conformation 
of  the  ERa  LBD  that  favors  coactivator  binding  and  one  of 
the  ERP  LBD  that  precludes  coactivator  binding.  However, 


because  it  lacks  a  bulky  side  chain  similar  to  those  of  OHT,  RAL 
and  ICI,  THC  must  antagonize  ERp  through  a  fundamentally 
different  mechanism  (Fig.  la).  Here  we  describe  structural  and 
additional  functional  characterization  of  the  effects  of  THC  on 
ERa  and  ERp.  These  studies  indicate  that  THC  exerts  its 
effects  on  the  AF-2  activity  of  the  two  ERs  by  stabilizing 
distinct  conformations  of  the  two  receptors  and  that  THC 
antagonizes  ERp  through  a  novel  mechanism  we  term  ‘passive 
antagonism’. 

Structures  of  the  THC-ER  LBD  complexes 

Crystals  of  ERa  LBD  bound  to  both  THC  and  a  peptide  contain¬ 
ing  the  second  NR  box  of  the  p  160  coactivator  GRIP1  (ref.  15), 
and  the  ERp  LBD  bound  to  THC  were  grown.  Both  structures 
were  determined  by  molecular  replacement.  The  ERa  complex 
structure  has  been  refined  to  an  R-factor  of  20.3% 
(Rfrcc  =  24.3%)  using  data  to  1.95  A,  and  the  ERp  complex  struc¬ 
ture  has  been  refined  to  an  R-factor  of  25.9%  (Rfrw  =  29.9%) 
using  data  to  2.95  A  (Table  1 ;  Figs  1  b,  3a). 

As  anticipated,  the  ERa  LBD  when  bound  to  THC  adopts  the 
same  conformation  as  it  does  when  bound  to  the  full  agonists  E2 
and  DES17’18  (Fig.  lb).  Helix  12  adopts  the  agonist-bound  con¬ 
formation,  which  promotes  NR  box  association.  In  contrast,  the 
binding  of  THC  to  ERp  LBD  does  not  stabilize  the  agonist- 
bound  conformation  of  helix  12  (Fig.  lb).  Surprisingly,  helix  12 
does  not  bind  to  the  region  of  the  coactivator  recognition  groove 
formed  by  helices  3,  4  and  5  as  it  does  in  the  OHT-ERa, 
RAL-ERa  and  RAL-ERp  complexes,  and  it  is  not  disordered  as 
it  is  in  the  ICI-ERp  complex  (Fig.  2a).  Instead,  helix  12  in  the 
THC-ERp  complex  adopts  an  orientation  most  similar  to  that 
observed  in  the  GEN-ERp  complex  (Fig.  2a).  In  both  complexes, 
helix  12  is  stabilized  in  this  conformation  by  hydrophobic  con¬ 
tacts  formed  by  Val  487,  Leu  491,  Met  494  and  Leu  495  with  the 
rest  of  the  LBD19  (data  not  shown). 
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Fig.  3  THC-!igand-binding  pocket  interactions, 
a,  Stereo  views  of  2F0-  Fc  electron  density  of  THC  in 
the  ERa  complex  (green)  and  in  the  ERp  complex 
(orange).  The  maps  were  calculated  to  1.95  A  and 
2.95  A  for  the  ERa  and  ERp  complexes,  respectively. 
Both  maps  were  contoured  at  1.0  a  and  calculated 
after  omitting  the  ligand,  b,  Schematic  diagrams  of 
the  interactions  between  THC  and  the  binding 
pockets  of  the  two  ERs.  Only  residues  or  waters  that 
form  hydrogen  bonds  with  THC  and/or  van  der 
Waals  contacts  with  THC  (4.2  A  cutoff)  are  depicted. 
Hydrogen  bonds  and  van  der  Waals  contacts  are 
represented  by  light  blue  and  gray  lines,  respective¬ 
ly.  Hydrogen  bond  distances  (A)  are  given.  Residues 
that  interact  only  with  THC  in  the  ERa  complex  are 
green,  and  those  that  interact  only  in  the  ERp  com¬ 
plex  are  blue-gray.  The  rings  and  the  individual 
atoms  of  THC  are  labeled,  c,  Stereo  view  of  the 
binding  pockets  of  the  two  ERs  bound  to  THC.  The 
structures  were  superimposed  using  LSQMAN34 
(1.1  A  r.m.s.  deviation  over  200  matched  Ca  atoms 
using  a  3.8  A  cutoff).  ERa  and  ERp  residues  are 
labeled  in  light  green  and  blue-gray,  respectively. 
Side  chain  atoms  are  colored  by  atom  type  (carbon 
(ERa)  =  light  green,  carbon  (ERp)  =  blue-gray,  nitro¬ 
gen  =  dark  blue,  oxygen  =  red  and  sulfur  =  yellow). 
THC  in  the  ERa  complex  is  colored  green,  and  THC  in 
the  ERp  complex  is  colored  orange.  Hydrogen  bonds 
are  depicted  as  dashed  black  bonds.  For  clarity, 
water  molecules  and  the  side  chains  of  Thr  347  and 
Leu  349  from  ERa,  and  Thr  299  and  Leu  301  from 
ERp  are  not  shown.  The  A  ring  and  A  ring  hydroxyl 
interact  comparably  with  the  equivalent  binding 
pocket  residues  from  the  two  ERs.  In  contrast,  the  A' 
ring  stabilizes  Met  421,  lie  424,  Gly  521,  His  524  and 
Leu  525  from  ERa  in  conformations  distinct  from 
those  of  lie  373,  lie  376,  Gly  472,  His  475  and  Leu  476 
from  ERp.  In  addition,  because  the  A'  ring  hydroxyl 
is  incorrectly  oriented,  it  fails  to  hydrogen  bond 
with  His  475  from  ERp  (as  it  does  with  the  imidazole 
of  His  524  from  ERa)  and,  instead,  weakly  hydrogen 
bonds  with  the  carbonyl  of  Gly  472. 


Ligands  and  helix  12  positioning 

Based  on  the  position  of  helix  12,  the  ERP  LBD  in  the  conforma¬ 
tion  stabilized  by  THC  and  GEN  should  be  incapable  of  inter¬ 
acting  with  coactivators  and,  hence,  transcriptionally  silent19. 
However,  GEN  and  THC  clearly  show  different  activities  on  ERp 
in  mammalian  cells:  GEN  is  a  partial  agonist9  and  THC  is  a  pure 
antagonist3’4.  How  can  the  structural  data  be  reconciled  with  the 
different  activities  of  these  compounds  in  transcriptional  assays? 

The  simplest  model  that  explains  these  data  is  based  on  two 
hypotheses.  First,  helix  12  in  the  unliganded  ERP  LBD  is  in  equi¬ 
librium  between  the  inactive  conformation  observed  in  the 
THC-  and  GEN-ERp  LBD  structures  and  the  active  agonist- 
bound  conformation.  Second,  rather  than  inducing  a  single 
static  conformation  of  helix  12,  ligands  affect  LBD  transcrip¬ 


tional  activity  by  shifting  this  equilibrium.  As  a  consequence, 
helix  12  in  the  unliganded  or  apo  ERp  LBD,  although  largely  in  a 
transcriptionally  silent  conformation,  may  also  adopt  a  tran¬ 
scriptionally  active  one.  Hence,  the  low  but  detectable  constitu¬ 
tive  AF- 2 -derived  transcriptional  activity  of  the  unliganded 
receptor9  results  from  the  modest  but  significant  ability  of  the 
apo  LBD  to  bind  coactivator.  Full  agonists,  such  as  E2  and  DES, 
shift  the  helix  12  conformational  equilibrium  in  favor  of  the 
active  conformation  and  stimulate  AF-2  activity  by  increasing 
the  affinity  of  the  LBD  for  coactivator.  Because  partial  agonists, 
such  as  GEN,  are  incapable  of  shifting  the  equilibrium  in  favor  of 
the  active  conformation  to  the  same  extent  as  full  agonists,  they 
are  able  only  to  incrementally  increase  the  affinity  of  the  receptor 
for  coactivator  and,  hence,  activate  AF-2  less  efficiently  than  full 
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agonists.  GEN  binding  should  permit  the  ER(3  LBD  to  sample 
both  inactive  and  active  conformations  of  helix  12  (possibly  to 
comparable  extents).  The  observation  of  the  inactive  conforma¬ 
tion  in  the  GEN  complex  structure  may  reflect  the  true  bias  of 
the  equilibrium  in  solution  or  it  may  result  from  the  influence  of 
the  crystallization  conditions  on  the  equilibrium. 

In  principle,  THC  could  antagonize  receptor  activity  by  shift¬ 
ing  the  equilibrium  to  favor  the  inactive  THC/GEN-bound  con¬ 
formation.  If  THC  binding  forces  the  receptor  to  sample  this 
conformation  more  than  it  does  in  the  absence  of  ligand,  the 
THC-receptor  complex  would  be  expected  to  have  an  even  lower 
affinity  for  coactivator  than  the  unliganded  receptor.  This  would 
agree  with  the  absence  of  transcriptional  activity  of  the 
THC-ERp  complex  and  the  conformation  of  helix  12  observed 
in  the  crystal. 

The  helix  12  conformational  equilibrium  would  be  difficult 
to  observe  directly.  Nonetheless,  the  predicted  actions  of  the 
various  ligands,  inferred  from  the  model  described  above,  sug¬ 
gest  a  testable  hypothesis.  The  binding  of  different  ligands 
should  modulate  the  affinity  of  the  ERp  LBD  for  coactivator 
such  that  E2-LBD  >  GEN-LBD  >  apo  LBD  >  THC-LBD.  The 
affinities  of  various  ligand-ERp  LBD  complexes  for  an 
Leu-X-X-Leu-Leu  motif-containing  peptide  were  directly  mea¬ 
sured  using  a  fluorescence  polarization-based  binding  assay21, 
and  the  expected  rank  order  of  affinities  was  observed  for  the 
different  complexes  (Fig.  2b).  Thus,  these  and  other  biochemi¬ 
cal  data5  suggest  that  helix  12  positioning  in  ERp,  as  in  the  case 
of  other  NRs33,  is  dictated  by  a  ligand-sensitive  conformational 
equilibrium. 

THC-ligand  binding  pocket  interactions 

When  ER  agonists  bind  to  ERa  (and  presumably  ERp),  they  favor 
the  active  conformation  of  helix  12  but  they  do  so  indirectly. 
These  ligands  form  hydrogen  bonds  and  van  der  Waals  contacts 
with  residues  from  helices  3,  6,  8  and  1 1  and  the  p-hairpin,  posi¬ 
tioning  these  residues  in  conformations  that  allow  their  inter¬ 
action  with  each  other  and  with  neighboring  residues 
surrounding  the  binding  pocket.  The  formation  of  these  cooper¬ 
ative  interactions  positions  residues  from  helices  3,  5  and  1 1  such 
that  they  create  a  predominantly  hydrophobic  binding  surface  for 
the  loop  connecting  helices  1 1  and  12  and  helix  12  itself7*18. 

THC  interacts  with  the  ERa  LBD  similar  to  the  full  agonists, 
DES  and  E2,  and  nucleates  all  of  the  polar  and  nonpolar  inter¬ 
actions  required  to  favor  the  positioning  of  helix  12  in  the  ago¬ 
nist-bound  conformation17*18  (Fig.  3 b,c).  For  example,  in  the 
ERa  complex,  the  A'  ring  forms  hydrogen  bonds  and/or  nonpo¬ 
lar  contacts  with  the  side  chains  of  His  524  and  Leu  525 
(Fig.  3b, c),  stabilizing  them  in  conformations  such  that  they 


Fig.  4  Ligand-binding  pocket  residues  influence  helix  12  positioning.  The 
structures  of  the  two  complexes  were  superimposed,  and  the  side  chains 
and  ligands  are  colored  as  in  Fig.  3c.  In  addition,  the  peptide  backbone 
of  ERa,  ERp  and  the  coactivator  peptide  are  colored  light  green,  blue- 
gray  and  purple,  respectively.  In  each  structure,  helices  11  and  12  are 
labeled  H1 1  and  H12,  and,  for  clarity,  the  ligands  have  been  made  semi¬ 
transparent.  The  side  chain  of  Met  479  from  ERp  is  disordered  and 
depicted  as  an  Ala  residue.  In  the  ERa  complex,  His  524,  Leu  525  and 
Met  528  are  oriented  by  the  THC  A'  ring  such  that  they  can  form  nonpo¬ 
lar  contacts  with  residues  from  helix  12  or  the  preceding  loop  (Val  533, 
Val  534  and  Leu  544).  Because  of  the  distinct  positioning  of  His  475  and 
Leu  476  (His  524  and  Leu  525  in  ERa)  stabilized  by  the  A'  ring,  the  equiv¬ 
alent  interactions  do  not  occur  in  the  ERp  complex. 


form  (either  directly  or  indirectly)  part  of  the  binding  site  for 
helix  12  or  the  preceding  loop.  The  side  chain  of  His  524  forms 
hydrophobic  contacts  with  the  side  chain  of  Met  528,  position¬ 
ing  it  such  that  it  can,  in  turn,  form  nonpolar  contacts  with  the 
side  chains  of  Val  533  and  Val  534  from  the  loop  between  helices 
11  and  12  (Fig.  4).  The  side  chain  of  Leu  525  packs  against  that  of 
Leu  544  from  helix  12  (Fig.  4). 

Although  it  binds  to  ERP  in  a  similar  overall  manner,  THC 
fails  to  stabilize  many  of  the  interactions  in  the  A'  ring  region  of 
the  binding  pocket  that  were  observed  in  the  ERa  complex 
(Fig.  3b, c).  In  the  ERp  complex,  the  B-ring  ethyl  group  adopts  a 
more  extended  conformation  than  in  the  ERa  complex,  causing 
the  tetrahydrochrysene  scaffold  to  be  tilted  by  5°  towards  the 
opening  formed  by  helices  3,  5  6  and  1 1,  and  away  from  the  floor 
of  the  binding  pocket.  Because  the  A'  ring  is  displaced  by  ~ 0.7  A 
from  its  position  in  the  ERa  complex  (Fig.  3c),  the  side  chains  of 
His  475  and  Leu  476  (which  form  only  nonpolar  contacts  with 
the  A'  ring)  are  positioned  ~1.6  and  2.3  A  away,  respectively, 
from  the  locations  of  their  counterparts  in  the  ERa  complex 
(His  524  and  Leu  525)  (Figs  3c,  4).  As  a  result,  the  side  chain  of 
His  475  fails  to  pack  against  that  of  Met  479  (equivalent  to 
Met  528  in  ERa),  causing  the  Met  479  side  chain  to  be  disor¬ 
dered.  The  backbone  atoms  of  Met  479  and  the  residues  flanking 
it  adopt  a  random  coil  conformation  (as  compared  to  the  helical 
conformation  adopted  by  their  counterparts  in  ERa),  and 
Met  479  is  located  2.4  A  distant  from  Met  528  (based  on  Ca 
atoms)  (Fig.  4).  Consequently,  both  Leu  476  and  Met  479  are  not 
positioned  appropriately  to  form  interactions  with  relevant 
residues  from  helix  12  and  the  preceding  loop  (Val  484,  Val  485 
and  Leu  495)  that  would  stabilize  the  active  conformation  of 
helix  12  (Fig.  4).  Therefore,  these  residues  destabilize  the  active 
conformation  of  helix  12  (Fig.  4).  Thus,  by  positioning  certain 
binding  pocket  residues  in  nonproductive  conformations,  the 
binding  of  THC  to  ER(3  disfavors  the  agonist-bound  conforma¬ 
tion  of  helix  12  and  shifts  the  equilibrium  towards  the  inactive 
conformation  of  helix  12. 

The  two  ligand-binding  pocket  residues  that  differ  between 
the  two  subtypes  may  explain  the  failure  of  THC  to  act  as  an  ERp 
agonist.  Met  336  and  lie  373  in  ERp  (equivalent  to  Leu  384  and 
Met  421  in  ERa)  are  positioned  on  opposite  faces  of  THC 
(Fig.  3b, c).  Modeling  (based  on  both  the  THC-ERa  and 
GEN-ERP  structures)  suggests  that  the  THC  A'  ring  would  have 
to  be  positioned  ~2.4  A  lower  in  the  ERP  binding  pocket  relative 
to  its  location  in  the  crystal  in  order  for  THC  to  both  interact 
with  Met  336  and  He  373  and  orient  His  475  and  Leu  476  in  pro¬ 
ductive  conformations  (data  not  shown).  This  lower  position 
should  be  highly  unfavorable  because  it  would  result  in  steric 
clashes  between  the  B'  ring  ethyl  group  and  the  side  chain  of 


362 


nature  structural  biology  •  volume  9  number  5  *  may  2002 


letters 


Table  1  Summary  of  Crystallographic  Statistics 


Data  collection 

Complex 

THC-ERa  LBD-GRIP1 

THC-ERP  LBD 

Space  group 

NR  Box  II  peptide 

P2i 

R3 

Resolution  (A)1 

1.95(2.02-1.95) 

2.95(3.06-2.95) 

Observations 

Total 

119,409 

64,540 

Unique 

34,533 

14,895 

Completeness  (%)’ 

97.9  (97.7) 

99.7  (100) 

Rsym  (%Y'2 

6.5(53.1) 

4.5  (30.2) 

<I>  /  <o  (I)>1 

18.8(2.9) 

23.3  (3.4) 

Refinement 

Number  of  atoms 
Protein 

3,830 

3,409 

Water 

166 

9 

Heterogen 

53 

48 

RcrystW1'3 

20.3  (29.6) 

25.9  (34.0) 

Rfree  (%)1*3 

24.3  (33.6) 

29.9(35.1) 

R.m.s.  deviation 

Bonds  (A) 

0.005 

0.011 

Angles  (°) 

1.08 

1.29 

Average  B -factor  (A2) 

38.5 

40.1 

’Values  in  parentheses  refer  to  the  highest  resolution  shell. 

2Rsym=  Ellh  -  <Ih>l  /  IIh,  where  Ih  is  the  integrated  intensity  of  a  given 
reflection  and  <Ih>  is  the  average  intensity  over  symmetry  equivalents. 
3Rcryst=  EIF0  -  Fci  / EIF0I,  where  F0  and  Fc  are  observed  and  calculated  ampli¬ 
tudes,  respectively.  Rfree  is  calculated  similarly  using  a  test  set  of  reflec¬ 
tions. 


Leu  298  at  floor  of  the  binding  pocket  (Fig.  3b, c).  THC  presum¬ 
ably  binds  in  the  alternative  mode  observed  in  the  crystal  to 
avoid  these  steric  clashes.  Consistent  with  this,  the  dimethyl  ana¬ 
log  of  THC,  which  should  be  less  sterically  hindered  from  adopt¬ 
ing  the  lower  position,  acts  as  a  weak  ERp  partial  agonist3. 

Antagonism  without  a  side  chain 

Because  the  positioning  of  the  side  chains  of  OHT,  RAL  and  ICI 
directly  or  ‘actively’  precludes  the  agonist-bound  conformation 
of  helix  12  by  steric  hindrance1344,  we  term  their  common  mech¬ 
anism  of  action  as  active  antagonism’.  Clearly,  THC  lacks  a 
bulky  side  chain,  and  in  its  complex  with  ERp,  helix  12  is  not 
sterically  precluded  from  adopting  the  agonist-bound  confor¬ 
mation  (Fig.  4).  Instead,  THC  antagonizes  ERp  by  stabilizing 
nonproductive  conformations  of  key  residues  in  the  ligand¬ 
binding  pocket,  thereby  disfavoring  the  equilibrium  to  the 
agonist-bound  conformation  of  helix  12  and  leading  to  stabiliza¬ 
tion  of  an  inactive  conformation  of  helix  12.  Thus,  we  call  the 
mechanism  of  antagonism  of  THC  ‘passive  antagonism’. 

Passive  antagonism  may  not  be  unique  to  THC  and  ERp. 
There  are  many  examples  of  NR  ligands  that  act  as  antagonists 
even  though  they  are  smaller  than  the  endogenous  agonists  of 
these  NRs.  The  synthetic  androgen  receptor  antagonist,  flu- 
tamide,  is  comparable  in  size  to  testosterone  and  does  not  pos¬ 
sess  an  obvious  moiety  that  would  act  as  an  antagonist  side 
chain22.  Similarly,  progesterone  is  smaller  than  aldosterone  but  is 
a  high  affinity  antagonist  of  the  mineralocorticoid  receptor23. 

Many  NRs  have  several  subtypes  that  possess  distinct  expres¬ 
sion  patterns  and  regulate  distinct  target  genes24.  Antagonists 
generated  through  the  addition  of  bulky  side  chains  to  agonist 
scaffolds  are  limited  to  being  antagonistic  on  one  or  more  sub¬ 


types  of  a  particular  NR1344.  The  passive  antagonism  mecha¬ 
nism,  as  revealed  here  through  direct  comparison  of  the  two 
THC-ER  LBD  complexes,  suggests  a  new  approach  to  achieving 
NR  antagonism.  Compounds  could  be  designed  to  selectively 
stabilize  the  inactive  conformations  of  certain  NR  subtypes  and 
the  active  conformations  of  others.  Such  ligands  may  exert  novel 
biological  and  therapeutic  effects. 

Methods 

Protein  expression  and  purification.  The  human  ERa  LBD 
(residues  297-554)  was  expressed,  carboxymethylated  and  purified 
as  described17.  For  crystallographic  studies,  the  human  ERp  LBD 
(residues  256-505)  was  expressed  as  an  N-terminally  His6-tagged 
fusion  protein  in  BL21(DE3)pLysS  cells  using  a  modified  pET-1 5b 
plasmid  (Novagen).  Bacterial  lysates  were  applied  to  an  estradiol- 
Sepharose  column,  and  the  bound  ERp  LBD  was  carboxymethylated 
with  20  mM  iodoacetic  acid.  Protein  was  eluted  with  30  pM  THC  in 
~50  ml  of  20  mM  Tris-HCI,  1  M  urea  and  10%  (v/v)  dimethyl  for- 
mamide  (DMF),  pH  8.1.  The  ERp  LBD  was  further  purified  by  ion 
exchange  chromatography  (Resource  Q,  Pharmacia).  Protein  sam¬ 
ples  were  analyzed  by  SDS-PAGE,  native  PAGE  and  electrospray  ion¬ 
ization  mass  spectrometry.  For  biochemical  studies,  residues 
214-530  of  human  ERp  were  expressed  as  a  fusion  with  glutathione- 
S-transferase  (GST)  in  BL21  cells  using  a  modified  pGEX-4T1  plasmid 
(Pharmacia).  The  GST  fusion  protein  was  bound  to  glutathione 
Sepharose  4  Fast  Flow  (Pharmacia)  and  eluted  with  glutathione  per 
the  manufacturer  instructions.  In  some  experiments,  the  protein 
was  further  purified  by  ion  exchange  chromatography  (HiTrap  Q, 
Pharmacia). 

Crystallization  and  data  collection.  Crystals  of  the  THC-ERa 
LBD-GRIP1  NR  box  II  peptide  complex  were  prepared  by  hanging 
drop  vapor  diffusion  at  19-21  °C.  Before  crystallization,  the 
THC-ERa  LBD  complex  was  incubated  overnight  with  a  four-fold 
molar  excess  of  the  GRIP1  NR  box  II  peptide.  Samples  (0.5  pi)  of  this 
solution  (5.0  mg  mM  protein)  were  mixed  with  3.5  pi  of  reservoir 
buffer  consisting  of  16%  (w/v)  PEG  4000,  53  mM  Tris-HCI,  pH  8.8, 
and  50  mM  MgCI2  and  suspended  over  wells  containing  800  pi  of 
the  reservoir  buffer.  The  crystals  lie  in  the  spacegroup  P21#  with  cell 
dimensions  a  =  54.55  A,  b  =  82.60  A,  c  =  59.04  A  and  p  =  1 1 1 .53°.  Two 
independent  THC-LBD-peptide  complexes  form  the  asymmetric 
unit.  A  crystal  was  transferred  to  a  cryo-solvent  solution  containing 
20%  (w/v)  PEG  4000,  15%  (v/v)  ethylene  glycol,  100  mM  Tris-HCI,  pH 
8.6,  and  100  mM  MgCI2  and  frozen  in  an  N2  stream  in  a  nylon  loop. 
Diffraction  data  were  measured  at  -170  °C  at  beamline  5.0.2  at  the 
Advanced  Light  Source  (ALS)  using  a  Quantum  4  CCD  camera  (Area 
Detector  Systems  Corp.)  at  a  wavelength  of  1.10  A.  Images  were 
processed  with  DENZO25,  and  the  integrated  intensities  were  scaled 
with  SCALEPACK25  using  the  default  -3  a  cutoff. 

Crystals  of  the  THC-ERp  LBD  complex  were  obtained  by  hanging 
drop  vapor  diffusion  at  21-23  °C.  Samples  (2  pi)  of  a  solution  of 
the  complex  (4.8  mg  mM)  were  combined  with  2  pi  samples  of  a 
reservoir  solution  containing  1.5-1.75  M  (NH4)2S04  and  100  mM 
sodium  acetate,  pH  4.8-5. 2,  and  suspended  over  wells  containing 
800  pi  of  reservoir  solution.  The  resulting  crystals  belong  to  the 
space  group  R3,  with  cell  parameters  a  =  b  =  99.14  A  and 
c  =  193.38  A  (hexagonal  setting).  The  asymmetric  unit  contains 
two  ERp  LBD  monomers  that  do  not  form  the  dimer  observed  in 
the  ERa  complex.  Instead,  each  of  the  two  LBDs  interacts  similarly 
with  symmetry-related  molecules  to  form  crystallographic  trimers. 
Before  data  collection,  a  single  crystal  was  transferred  to  a  stabi¬ 
lizing  solution  (1.8  M  (NH4)2S04,  100  mM  NaCI,  100  mM  sodium 
acetate,  pH  4.5,  and  10  pM  THC).  The  crystal  was  then  sequentially 
transferred  at  30  min  intervals  through  a  series  of  solutions  con¬ 
sisting  of  the  stabilizing  solution  supplemented  with  increasing 
concentrations  of  ethylene  glycol  (1%  (v/v)  increments)  to  a  final 
concentration  of  15%  (v/v).  The  crystal  was  then  flash  frozen  in  an 
N2  stream  in  a  nylon  loop.  Diffraction  data  were  measured  at 
-170  °C  at  beamline  5.0.2  at  the  ALS  using  a  Quantum  4  CCD 
camera  at  a  wavelength  of  1.07  A.  Images  were  processed 
with  DENZO,  and  the  integrated  intensities  were  scaled  with 
SCALEPACK  using  the  default  -3  a  cutoff. 
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Structure  determination  and  refinement.  The  model  of  the 
DES-ERa  LBD-GRIP1  NR  box  II  peptide  complex  (3ERD)  was  modi¬ 
fied  by  truncating  the  side  chains  of  the  ligand-binding  pocket  of 
both  monomers  to  Ala  and  by  removing  all  ligands,  waters,  car- 
boxymethyl  groups  and  ions.  After  subjecting  this  edited  model  to 
rigid  body  refinement  in  REFMAC26  the  missing  parts  of  the  model 
were  built.  The  rest  of  the  model  was  corrected  using  MOLOC27  and 
two-fold  averaged  maps  generated  with  DM26.  All  masks  for  aver¬ 
aging  were  generated  using  MAMA28,  and  PHASES29  and  the  CCP4 
suite26  were  used  for  the  generation  of  structure  factors  and  the  cal¬ 
culation  of  weights.  Initially,  positional  refinement  was  performed 
using  REFMAC.  At  later  stages,  the  model  was  refined  using  the  sim¬ 
ulated  annealing,  positional  and  B-factor  refinement  protocols  in 
CNS30.  All  B-factors  were  refined  isotropically.  Anisotropic  scaling 
and  a  bulk  solvent  correction  were  used  throughout  refinement. 
The  R uee  set  contains  a  random  sample  of  6.5%  of  all  data.  All  data 
between  47  and  1.95  A  (with  no  a  cutoff)  were  used.  The  current 
model  is  composed  of  residues  305-459  and  469-547  of  monomer 
A,  residues  305-460  and  472-547  of  monomer  B,  residues  687-696 
of  peptide  A,  residues  686-695  of  peptide  B,  two  ligand  molecules, 
166  water  molecules,  one  carboxymethyl  group  and  one  chloride 
ion.  According  to  PROCHECK26,  96.6%  of  all  residues  in  the  model 
are  in  the  core  regions  of  the  Ramachandran  plot,  and  none  are  in 
the  disallowed  regions. 

The  intensities  within  the  ERp  data  set  fall  off  very  rapidly  as  a 
function  of  resolution  (Wilson  B-factor  -90  A2  calculated  from 
4.5  to  2.95  A).  To  increase  the  contribution  of  the  higher  resolution 
terms,  the  data  were  sharpened  with  a  correction  factor  of  -55  A2. 
Maps  calculated  using  this  sharpened  data  revealed  higher  resolu¬ 
tion  features,  such  as  improved  side  chain  density31.  All  subsequent 
manipulations  were  performed  using  this  sharpened  data. 

The  two  LBDs  in  the  asymmetric  unit  were  located  by  molecular 
replacement  in  AmoRe26  and  TFFC26.  The  search  model  was  con¬ 
structed  by  overlapping  the  models  of  five  ERa  LBD  complexes  (PDB 
entries  1A52,  1ERE,  1ERR,  3ERD  and  3ERT)  and  setting  the  occupan¬ 
cies  of  each  model  to  20%  (R-factor  =  55.3%  and  correlation  coeffi¬ 
cient  =  55.9%  after  placement  of  both  monomers).  The  model  was 
then  built  using  MOLOC  and  two-fold  averaged  maps  generated 
with  DM.  MAMA  was  used  for  all  mask  manipulations,  and  PHASES 
and  the  CCP4  suite  were  used  for  the  generation  of  structure  factors 
and  the  calculation  of  weights.  Refinement  was  performed  initially 
using  the  positional  refinement  protocols  in  REFMAC  and,  later, 
using  the  simulated  annealing,  positional  and  B-factor  refinement 
protocols  in  CNS.  All  B-factors  were  refined  isotropically. 
Anisotropic  scaling,  a  bulk  solvent  correction,  and  tight  noncrystal- 
lographic  symmetry  restraints  were  used  throughout  refinement. 
The  R,ree  set  contains  a  random  sample  of  5%  of  all  data.  All  data 
between  49  and  2.95  A  (with  no  a  cutoff)  were  used.  The  current 
model  is  composed  of  residues  261-284,  290-408  and  413-501  of 
monomer  A,  the  last  residue  of  the  affinity  tag  and  residues 
256-  280,  294-410,  413-479  and  485-497  of  monomer  B,  two  ligand 
molecules,  and  nine  water  molecules.  According  to  PROCHECK26, 
92.5%  of  all  residues  in  the  model  are  in  the  core  regions  of  the 
Ramachandran  plot,  and  none  are  in  the  disallowed  regions. 

Peptide  binding  assay.  Samples  of  the  GST-ERp  LBD  fusion  pro¬ 
tein,  in  the  absence  or  presence  of  the  various  ligands  (at  saturating 
concentrations),  were  mixed  with  an  Leu-X-X-Leu-Leu  motif-con¬ 
taining  peptide  (lle-Leu-Arg-Lys-Leu-Leu-GIn-Glu),  which  had  been 
N-terminally  labeled  with  rhodamine  (reaction  buffer  =  10  mM 
HEPES-Na,  150  mM  NaCI,  2  mM  MgCI2,  1  mM  EDTA  and  100  pg  mM 
BSA,  pH  7.9,  and  [peptide]  =  5  nM).  The  binding  reactions  were 
incubated  at  ambient  temperature  with  shaking  for  1  h  in  96-well 
plates  (Whatman),  and  the  fluorescence  polarization  was  measured 
on  an  Analyst  reader  (Molecular  Devices).  Kd  values  were  generated 


from  triplicate  assays  by  nonlinear  least  squares  analysis  (Prism, 
GraphPad  Software). 

Coordinates.  Coordinates  of  the  structures  have  been  deposited  in 
the  Protein  Data  Bank  (accession  codes  1L2I  and  1L2J  for  the  ERa 
and  ERp  complexes,  respectively). 
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